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Abstract 

We study the magnetization distribution of the Ising model on two regular 
fractals (a hierarchical lattice, the regular simplex) and percolation clusters at the 
percolation threshold in a two dimensional imbedding space. In all these cases, 
the only fixed point is T = 0. In the case of the two regular fractals, we show 
that the magnetization distribution is non trivial below T* ~ A*/n, with n the 
number of iterations, and A* related to the order of ramification. The cross-over 
temperature T* is to be compared with the glass cross-over temperature T g ~ A g /n. 
An estimation of the ratio T*/T g yields an estimation of the order of ramification 
of bidimensional percolation clusters at the threshold (C = 2.3 ± 0.2). 



1 Introduction 



We consider the problem of Ising spin systems on a family of fractals with a zero temperature 
fixed point. More especially, we consider the case of a hierarchical lattice, the regular simplex 
(regular fractals), and bidimensional percolation clusters at the percolation threshold. This 
problem has its origin in the pioner work of Mandelbrot [|TJ. The problem is also of an 
experimental interest because of the existence of random magnets with a fractal structure (see 
for instance M). As a byproduct, we derive a measure of the ramification order of percolating 
clusters. To do so, we compare to cross-over temperatures: 1) the glass cross-over temperature 
T g 2) the magnetization distribution cross-over temperature T*. The dominant behavior of 
the glass cross-over temperature was discovered in the 80's [[J [f|. This quantity is related 
to the correlation length of the magnet. Below T g , the correlation length is larger than the 
system size, leading to a linear regime of the logarithm of the relaxation time as a function of 
the inverse temperature. Above T g , the correlation length is smaller than the linear size of the 
cluster, leading to a parabolic behavior for the logarithm of the relaxation time as a function 
of the inverse temperature. This behavior is well established on a numerical basis ||. The 
second cross-over temperature T* was not much studied in the past. Below T* [which is inverse 
proportional to the logarithm of the number of sites], the magnetization distribution is non 
Gaussian. These non Gaussian magnetization distributions also accur on the ferromagnetic 
Cayley tree ||, however with a different decay with the system size. More surprising is the 
appearance of similar magnetization distributions in magnetized spin glasses 0. However, in 
real spin glasses, this non Gaussian magnetization distribution is not expected to vanish in 
the thermodynamic limit as it is the case on ferromagnetic systems. We first study the case of 
regular fractals (a hierarchical lattice and the regular simplex). We use the Swendsen algorithm 
to generate equilibrium states ||, and thus to calculate the magnetization distributions. In 
the case of the regular simplex, it is possible to determine the geometric origin of the different 
local maxima of the magnetization distribution at low temperatures. In section we make a 
general argument to estimate the magnetization distribution cross-over scale, which is shown 
to be inverse proportional to the number of sites, and proportional to the ramification order 
C . In the case of percolating clusters, we calculate the magnetization distribution at low 
temperature. Averaging over the geometry, we can see the existence of finite size effects: the 
spin system is more ferromagnetic if one lowers the number of generations. We also studied 
the correlation between different equilibrium states, which gives us an other way to estimate 
T*. Finally, we calculate the ration T g /T*, which yields an estimation of the ramification 
order of bidimensional percolation clusters. We end up with some concluding remarks. 
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2 A Hierarchical Lattice 



The interest of hierarchical lattices is that the thermodynamics of the Ising model is exactly 
soluble on these structures. We consider the case of the hierarchical lattice built with the rules 
of figure |l|. The coordination of every site is bounded above, which is not always the case for 
hierarchical lattices (for instance for the hierarchical diamond, in which case the coordination 
diverges in the thermodynamic limit). The number of sites of the hierarchical lattice of figure 
fj] is N n = 4.6 n /5 + 6/5, whereas the end to end distance is L n = 4 n_1 , leading to the fractal 
dimension d = ln6/ ln4. 

An interesting aspect of hierarchical lattices is that their Yang and Lee zero can be cal- 



culated easily in the temperature plane. This was done in |T0[ in the case of the hierarchical 
diamond. The Yang and Lee zeros set is a Julia set. For comparison with the results of [|T0[ |, 
we also calculate the Yang and Lee zeros of our hierarchical lattice. First, we need to derive 
recursion relations for the partition function. 

2.1 Recursion for the Partition Function 

The technique (standard for hierarchical lattices) consists in carrying out the trace over the 
the spin variables at the smallest scales, and to derive the renormalization of the temperature 
and magnetic field. The partition function of the hierarchical lattice with one generation is 

Z 1 (Y i £') = ^ e / 3J ( (J 0S+(^0+(73)((Ti+ ( 72))+(T32' e /3//(tT0 + ...+CT3) ^ 

CT0,<Tl,<T2,O'3 

where the spins are pictured on figure [1|. A straightforward calculation leads to 

Zi (£,£') = 2 (e 4/3J cosh (/3 J(£ + £') + 4/3#) + e" 4/3J cosh/3J(£ + £') 

+ cosh (p J(E - £') + 2/3H) + cosh (/3J(£ - £') - 2/3H) (2) 
+2 cosh/3 J(£ - £') + 2 cosh (/3J(£ + £') + 2(3H)) . 
The partition function -Zi(£, £') can be rewritten under the form 

Z 1 (£,£ / )=Are /3jES 'e /3 ^ E+E '). (3) 

We have introduced here three parameters A/", J and H for three distinct partition functions 
Z(+, +), Z(—,—) and Z(+,—) = Z(—,+). The partition function can thus be consistently 
be brought under the form (|3]). The identification of (0) and (|3]) leads to 

Af = fz 1 (+,+)Z 1 (+, -)%(-,-)) 1/4 (4) 



e*pW = f£±di. (6) 
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In the case of a zero magnetic field, we obtain 



with 



and 



Z n (J)=Af* n - 1 Z n _ 1 (j), (7) 

M 2 = 8 (cosh 6(3 J + 3 cosh 2(3 J + 4) (cosh A(3 J + 2 cosh 2(3 J + 1) , (8) 

2f3 j cosh 6/5 J + 3 cosh 2(3 J + 4 
6 = 2 (cosh 4(3 J + 2 cosh 2(3 J + 1) ' ^ 

The renormalization equation (|9|) has only one fixed point: J = 0, so that the spin system is 
always disordered at any finite temperature in the thermodynamic limit. 

2.2 Digression: Yang and Lee zeros in the Temperature Plane 

We now calculate the Yang and Lee zeros in the temperature plane. This is not the main 
subject of this paper, but it is interesting to compare with existing results on the hierarchical 



diamond, where the Yang and Lee zero in the temperature plane form a Julia set |10| . We 
work in the plane of the variable z = exp2(3J and invert the relation @, which leads to a 
polynomial of degree 6: 

z 6 - 2~zz h + (3 - Az)z^ + 4(2 - z)z 3 + (3 - 4z)z 2 - 2zx + 1 = 0, (10) 



If z is a solution of ( 10]) , 1/z is a solution too, so that the set of zeros of the partition function 
is invariant under the inversion z — > 1/z. In order to compute the set of zeros of the partition 
function, we use the methods of [[Uj], that is we start from a given point in the complex plane, 
we calculate the zeros of (|9]), choose one solution among the six zeros at random, and reiterate. 
The zeros are computed using a Laguerre's method [IT]. The resulting set is independent on 



the initial value of z, as far as one eliminates the first zeros. The result is plotted on figure 
0. The set of zeros are concentrated on lines, which intersect the real axis in the vicinity of 
the point z — 0. A zoom in the vicinity of z = reveals that the density of zeros vanishes 
as one approaches the origin (see figure §). Contrary to the case of the hierarchical diamond, 
the Lee and Yang zeros are concentrated on lines in the z = exp 2(3 J plane, whereas the Julia 
set of the hierarchical diamond is a fractal. 

2.3 Glass Cross-over Temperature 

The glass temperature cross-over T g is related to the decay of the correlation function from one 
end of the fractal to the other || []5J ||. Since the Ising model is soluble on the hierarchical 
lattice, it is easy to get the behavior of the correlation length between the extremal sites as 



3 



a function of the number n of iterations and the temperature. We call x the dimensionless 
inverse temperature x = (3 J. In the low temperature regime, we get from equation @ 
x n+ i ~ x n — In 2/2, where x n is the value of x after n renormalization steps. In the low 
temperature regime, x n+ \ ~ 1x A n . The limit between the high and low temperature regimes 
corresponds to the minimum k of 2x 4 — a; + In 2/2. We find K ~ 0.59. The correlation between 
the two extremal sites is nothing but 

where 

Z\p\ oc J^eF*, (12) 

so that the correlation is simply (EE') = tanh/5J n , which decreases from unity to zero as n 
increases from 1 to +oo. The cross-over temperature T g corresponds to x n > k and x n+ i < k, 
leading to the dominant behavior of the glass cross-over temperature 

T » = 5^2- < 13 > 

2.4 Magnetization Distribution 

The magnetization distribution becomes non gaussian below a temperature T* . This is due to 
the existence of macroscopic domains that are weakly connected ot the rest of the structure. 
The temperature cross-over T* is evaluated in section f|in a more general context, so that we 
do not reproduce the argument here in the case of the hierarchical lattice. The result in the 
case of the hierarchical lattice is 

4 7 

T* = -p-:, (14) 
nmo 

where we have kept only the dominant behavior for large n. The results for the magnetization 
distribution are plotted on figure § for n — 4. Using (0) for n = 4, we get T* ~ 0.55, which 
is consistent with the data of figure f|. Notice that, in this case, T g < T*. 

3 Regular Simplex 

The regular simplex recursion is plotted on figure [5| The number of sites if N n = 3 n , and the 
corresponding distance is L n = 2 n , so that the fractal dimension is d = In 3/ In 2. 



4 



3.1 Recursion for the Partition Function 

The partition function in a zero external field can be calculated recursively by calculating the 
trace over the spins at the deepest generation. With the notations of figure [5], the partition 
function is 

Z(Ei,E 2 ,E 3 ) = exp (/3J(Si(o-i2 + 0-13) + S 2 (o- 2 i + 0-23 + S 3 (o- 3 i + cr 32 )) (15) 
exp (/3J(cr 3 icr 3 2 + cr 2 iO"23 + C12C13)). (16) 

We look for Z under the form 

Z(E X , £ 2 , S 3 ) = AT exp (p J(E X E 2 + £ 2 £ 3 + £^3)) • (17) 

The renormalized form of Z has two parameters (N and J), for two distinct equations: one 
for (Ei, E 2 , S 3 ) G <Si, and the other for (£ x , E 2 , E 3 ) G S 2 , with 

51 = {(+,+,+),(-,--)} (18) 

5 2 = {(+,+,-),(+,-+),(-+,+),(--+),(-+,-),(+,--)}■ (19) 

The parameters TV and J are simply given by 

AT = (Z(+,+,+)Z(+,+,-)) 1/4 (20) 

Carrying out the trace over the a variables leads to 

Z(+,+,+) = 27e" 3/3J + 27e /3J + 9e 5/3J + e 9/3J (22) 
Z(+,+,-) = 3e~ 7(3J + 19e~ 3/3J + 33e pj + 9e 5(3J . (23) 

In the high temperature regime, the coupling constant renormalization is x ~ x 3 and, in the 
low temperature regime, x ~ x — (ln3)/2, where x = (3 J. The cross-over between the high 
and low temperature regime is x = K, with k the minimum of x 3 — x + (ln3)/2 (/t ~ 0.58). 
The correlation function is 

^>=mk^ =% (24) 

with 

Z\ji] oc ex P (^n(SiE 2 + E!E 3 + E 2 E 3 ) exp (//E 1 E 2 ) = — — — — , (25) 

from what we deduce the dominant behavior of the cross-over temperature scale 

T 9 = ( 26 ) 
nm3 
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3.2 Magnetization Distribution 



The magnetization distribution is plotted on figure |6] for various temperatures. Again, we do 
not reproduce the calculation of the magnetization cross-over temperature T* since a general 
argument in developed in section |j. The result for the dominant behavior of T* in the limit 
of large n is 

6 7 

T* = 4-r. (27) 



nln3 

In the case n = 7 (which corresponds to figure ^), we have T* ~ 0.8, which is consistent with 
the data of figure Notice also that, in the low temperature regime, when the non gaussian 
structure of the magnetization distribution is fully developped, it is possible to identify the 
origin of the different local maxima of P{m). To do so, we carry out a low temperature 
expansion of the magnetization distribution around the broken symmetry state m — 1, up to 
the order /z 3 , where n = e _/3J /(e /3J + e _/3J ) ~ e~ 2/3J . The result i 



P(m) 



is 



n-l 

1 - 72/i 2 - A^)5{\ - m) + 3/i 2 J2 s ( m - 1 + 2 - 3P ~" 

n-2 

+// ^(3 n ~ p - 3)5(m - 1 + 2.3 p " n ) + 6 £ 5(1 



2M X 



(28) 



+ o(/x 3 ). (29) 



t4 is a normalization constant and the set A is 

A = {M G (1, [JV/2]), 3A;, 3/, < I < k, M = 3 k ± 3'}. 



(30) 



We can recognize on figure |6] the contributions from m 
1/3 ± 1/27, 1/3 ± 1/9, 1/9 ± 1/27. 



1/3, 1/9, 1/27 but also from m 



4 General Argument for Regular Fractals with a Zero 
Temperature Transition 

We first recall what happens on the Cayley tree 0, where there exists also a glassy cross-over, 
with T g (x J I Inn, where n is the number of generations of the tree, and also a crossover in 
the magnetization distribution, with T* oc J / In n. It is possible to derive exact recursion 
relations for the average magnetization, to solve the recurence, and to obtain the temperature 
cross-over T*. The cross-over temperature is interpreted as follows: below T*, there exists less 
than one broken bound on a path connecting the center of the tree to the leaves, which leads 
to the criterium for the magnetization distribution cross-over temperature: 

n/i ~ 1. (31) 
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Taking the logarithm leads to T* oc J / In n. 

In the case of fractals, there are not bifurcations, so that the criterium corresponding to 



(HD is 

/i C iV n ~ 1, (32) 

which means that, at the cross-over, there exists only C broken bounds on the graph, where 
C is the number of bounds that need to be cut to isolate a p generations fractals in the bulk 
of an generations fractal. This is the ramification order. We thus obtain 

r - < 33 > 

In order to check this equation, we come back to the case of the hierarchical lattice, where we 
obtained recursively the partition function in an external magnetic field. We impose a small 
magnetic field, and calculate the partition function in the presence of the small magnetic field, 
and deduce the average magnetization numerically as a function of the number of generations. 
For a given number of generation, we look for the temperature for which the magnetization 
is a given constant. We thus obtain the curve /3*(n), which is expected to be a straight line 
from equation (|33|). This is indeed the shown on figure |7], which validates the previous 

heuristic argument. 

5 Percolation Clusters at the Threshold 

5.1 Glass cross-over temperature 

Percolation clusters at the percolation threshold were studied in the past as an exemple of criti- 
cal dynamics, with a glassy-like dynamics below T g @ ■ The glass cross-over temperature 
T g is evaluated in || with the argument that below T g , the correlation length [|EJ 

£t = exp (^) (34) 
is inferior to the linear size of the cluster. One gets easily || 

T - (35) 
where N is the number of sites of the percolation cluster. 

5.2 Magnetization Distribution 

The magnetization distribution of a single cluster at low temperatures looks like the one of 
regular fractals, apart from the fact that the localization of the local minima depends on 
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the geometry and is not easy to localize. For instance, we generated a percolation cluster at 
threshold in a bidimensional imbeding space (p c = 0.593). The corresponding low temper- 
ature magnetization distribution is plotted on figure |] for various temperatures. We show 
on figure |9] two configurations belonging to the local maximum with a normalized number of 
up spins N^/N ~ 0.38. The interpretation of the magnetization distribution thus depends 
crucially on the geometry. It is therefore legitimate to ask whether the local maxima structure 



persists after averaging over the geometry. The answer is 'no': we plotted on figure |10| the 
magnetization distribution averaged over the geometry for different temperature. We observe 
no local maxima. If the temperature increases, the weight of the ferromagnetic maximum 
decreases, whereas the tail of the distribution increases, which is consistent with the fact that 
the magnetization distribution is gaussian in the high temperature regime. 

Moreover, by averaging over the geometry, we can adress the question of finite size effects. 
We plotted on figure [II] the distribution of magnetization for two different sizes at the same 
temperature. We see that the small clusters have a more pronounced ferromagnetic peak, 
which is in agreement with the fact that T* decreases with the system size. 

5.3 Correlations 

We wish to analyze the correlation between the different spin configuration at a given tem- 
perature. To do so, we introduce n replica of the spin system. We call q\ a ^ G {0, 1} the value 
of the Potts variables in the replica a (1 < a < n), and we consider the following correlation: 

1 N 9 

R N (q, T) = - E 7 77 E fffV (36) 
qf^n(n- 1) ^ 

for a given geometry with a given magnetization (all the replica have the same magnetization). 
The variable q is the average magnetization: 

1 N 

v«,g = ^E# } - (37) 

i=l 

The correlation Rn{q, T) in a sector of given magnetization is nothing but a kind of Edwards- 
Anderson order parameter in a subspace of given magnetization. In the high temperature 
regime, qf is with a probability 1/2 and 1 with a probability 1/2, so that the average of 
qf^q^ is 1/4, so that 

R N (q,t) = ^- (38) 
if T > T*. We are interested in the variations of R for a given magnetization as a function 



of the temperature. These variations are plotted on figure O for a given geometry. The 
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magnetization is q = N^/N ~ 0.38. The four different curves are plotted for neighboring 
magnetizations. As expected, R tends to 0.38 above T*. The estimation of T* from this 
technique is coherent with the estimation from the analysis of the magnetization distribution 
(see figure ||). 

5.4 Evaluation of the Ramification Order of Percolating Clusters 

We know from two different methods that the cross-over temperature T* for N = 2660 sites is 
about T* ~ 0.6. From this data, we can deduce the ramification order C, which, in the case 
of regular fractals was interpreted as the number of links to be cut to isolate a p generations 
fractal inside a n generations fractal (n > p). We deduce from (|33|) and fl35|) that 

T - = (39) 

On the other hand, we know that v v ~ 1.3 in a bidimensional imbedding space [jnj, and that 
d = 1.78 ± 0.02 0. With N = 2660 sites, we evaluate T g ~ 0.6 from flU). Using (||), we 
get C = 2.3 ± 0.2. The incertainty comes from the incertainty in the location of the cross-over 
temperature T*. This result is in agreement with the fact that the order of ramification of 
percolating clusters is believed to be finite and superior to two (percolating clusters are not 
quasi one dimensional) [pi 



6 Conclusion 

We have studied the magnetization distribution on fractals with a zero temperature transition. 
Below a temperature inverse proportional to the logarithm of the number of sites, and propor- 
tional to the order of ramification, the magnetization is not trivial, with local maxima. This 
is due to the fact that is is possible to cut the graph into many large parts, by breaking only 
a small number of links. We have introduced a magnetization-dependent Edwards- Anderson 
order parameter, the variations of which are correlated with the magnetization distribution. 
It would be interesting to have an analytical proof of this fact. We have shown in that the 
cross-over temperature T* of percolating clusters is related to the order of ramification of the 
structure. We find C = 2.3 ± 0.2 for bidimensional percolation clusters. The cross-over tem- 
perature T* can be estimated in two different ways: 1) by calculating directly the probability 
distribution at different temperatures 2) by calculating the correlation between equilibrium 
states with a given magnetization. Our result is in agreement with the fact that the percola- 
tion cluster at threshold is not quasi one dimensional (which was known a long time ago !). 
However, it is clear that our method to determine the order of ramification is not accurate, 
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since one has to estimate a cross-over temperature, which is determined up to a certain incer- 
tainty. It would be interesting to generalize this study to magnetized spin glass phases. This 
will be done in a near future. 

I acknowledge J.C. Angles d'Auriac who lent me his Swendsen program, and I thank B. 
Dougot for comments on the manuscript. 
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Figure captions 



Figure 1: 

Construction of the hierarchical lattice. At the n = step, we start with two sites connected. 
After one step, the lattice contains 6 sites. At each step, the links are transforms as in the 
n = — ■> n = 1 transformation. 



Figure 2: 

Set of zeros of the partition function in the plane z = exp2(3J We have plotted 39000 zeros 
calculated by the procedure described in the text. 1000 zeros have been eliminated. 



Figure 3: 

Zoom of the set of zeros in the vicinity of the origin, for 19000 zeros of the partition function. 
The density of points corresponds to the density of zeros. The density of zeros vanishes in the 
vicinity of the origin. 



Figure 4: 

Magnetization distribution for the hierarchical lattice with 4 generations. For clarity, the plots 
have been shifted along the y axis. 750000 iterations of Swendsen algorithm were carried out. 



Figure 5: 

Construction of the regular simplex. 



Figure 6: 

Magnetization distribution for the regular simplex with 7 generations. For clarity, the plots 
have been shifted along the y axis. 750000 iterations of Swendsen algorithm were carried out. 



Figure 7: 

Inverse temperature cross-over as a function of the number of generations for the hierarchical 
lattice. The cross-over is such that the average magnetization is some constant, the two 
extremal spin being frozen in the up direction. The constant is taken equal to 0.02. 
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Figure 8: 

Magnetization distribution for a percolation cluster of size iV = 2660 sites (this cluster is 
pictured on figure |.The temperature is T = 0.34,0.42,0.50,0.58,0.66,0.74. 750000 iteration 
of the Lanczos algorithm were carried out. 



Figure 9: 

Two spin configurations in the peak N-^/N ~ 0.38. 



Figure 10: 

Magnetization distribution averaged over the geometry. The clusters are generated in a 
60x60 box, and 100000 steps of Swendsen algorithm are carried out. For T = 0.4(0.3,0.5), 
106(200, 58) different geometries were generated. Notice the semi-log scale. 



Figure 11: 

Magnetization distribution averaged over the geometry. The clusters are generated in a 60x60 
box and a 25x25 box, and 100000 steps of Swendsen algorithm are carried out. The tempera- 
ture is T = 0.25. 



Figure 12: 

Variations of the correlation R as a function of the temperature for the cluster of figure 
The average number of up spins is chosen to be around 0.38. 
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